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Preface 


This  report  describes  the  soil-to-structure-to-limestone  foundation  interaction 
analysis  used  to  assess  potential  wall  performance  of  the  new  roller-compacted 
concrete  north  lock  wall  at  McAlpine  Locks.  Funding  for  this  study  was 
provided  by  the  U.S.  Army  Engineer  District  (USAED),  Louisville.  Mr.  Larry 
Dalton  and  Ms.  Monica  Greenwell,  Design  Branch,  USAED,  Louisville,  were 
the  project  monitors  under  the  general  supervision  of  Mr.  Byron  McClellan, 
Chief,  Design  Branch. 

This  work  was  performed  at  the  U.S.  Army  Engineer  Waterways  Experiment 
Station  (WES)  by  Dr.  Robert  M.  Ebeling,  Computer-Aided  Engineering  Division 
(CAED),  Information  Technology  Laboratory  (ITL),  and  by  Mr.  Ronald  E. 

Wahl,  Soil  and  Rock  Mechanics  Division,  Soil  Research  Center,  Geotechnical 
Laboratory.  Mr.  John  Hendricks,  ITL,  provided  invaluable  assistance  in  process¬ 
ing  the  results  of  the  computer  analyses  and  preparing  the  figures  for  this  report. 
The  authors  of  this  report  are  also  indebted  to  Mr.  Kevin  Abraham,  ITL,  for 
rendering  his  timely  assistance  in  the  development  of  the  finite  element  mesh  of 
the  lock.  This  report  was  prepared  by  Dr.  Ebeling  and  Mr.  Wahl.  Review 
commentary  was  also  provided  by  Mr.  Dalton  and  Ms.  Greenwell.  The  work 
was  accomplished  under  the  general  direction  of  Mr.  H.  Wayne  Jones,  Chief, 
CAED,  ITL,  and  Dr.  N.  Radhakrishnan,  Director,  ITL. 

At  the  time  of  the  publication  of  this  report,  Director  of  WES  was 
Dr.  Robert  W.  Whalin.  Commander  was  COL  Bruce  K.  Howard,  EN. 


The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or 
promotional  purposes.  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 


Soil-Structure-Foundation 
Interaction  Analysis  of  New 
Roller-Compacted  Concrete 
North  Lock  Wall  at  McAlpine 
Locks 


Background 

A  new  lock  is  to  be  constructed  at  the  McAlpine  Locks  on  the  Ohio  River  at 
Louisville,  KY.  The  new  lock  chamber  will  be  defined  by  rock-founded,  gravity, 
earth-retaining  monoliths.  The  design  of  these  new  lock  walls  was  completed  by 
Louisville  and  Portland  District  personnel.  One  of  the  alternatives  is  to  construct 
the  walls  using  roller-compacted  concrete  (RCC)  with  a  concrete  facing,  with  the 
culverts  located  in  the  floor  of  the  lock  chamber.  The  objective  of  this  study  was  to 
assess  the  soil-to-structure-to-foundation  interactions  of  a  typical  RCC  lock  wall 
two-dimensional  (2-D)  section. 

This  report  describes  the  results  of  a  complete  soil-structure  interaction  (SSI) 
analysis  of  a  typical  section  of  the  north  wall  of  the  new  McAlpine  Lock  monolith 
using  the  backfill  placement  method  of  analysis  incorporated  in  the  finite  element 
computer  program  SOILSTRUCT- ALPHA.  SOILSTRUCT  is  used  to  understand 
the  sometimes  complex  interactions  among  the  lock,  the  backfill,  and  the  founda¬ 
tion  rock  strata.  The  ALPHA  version  of  SOILSTRUCT  (Ebeling,  Duncan,  and 
Clough  1990)  has  been  developed  to  analyze  gravity  lock  walls  like  the  McAlpine 
lock  walls. 
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Section  Analyzed 

The  north  wall  of  the  new  RCC  McAlpine  Lock  is  shown  in  Figure  1 .  The  lock 
wall  being  analyzed  is  73  ft  (22.25  m)  high  (above  elevation  3701 ),  47  ft  (14.33  m) 
wide  and  retains  55  ft  (21  m)  of  well-compacted,  dense  granular  back-fill.  The  base 
of  the  new  lock  (at  el  370)  is  9.5  ft  (2.9  m)  above  the  floor  of  the  new  lock  chamber. 
The  postconstruction  water  table  in  the  backfill  is  assigned  el  395,  and  the  pool 
elevation  is  at  el  383  in  the  complete  SSI  analysis.  The  section  analyzed  is 
representative  of  the  north  wall  of  the  new  lock  at  station  25+00.  The  complete  2-D 
section  modeled  in  the  analysis  is  shown  in  Figure  2.  The  section  includes  two  of 
the  existing  lock  walls,  the  existing  backfill  and  shale  foundation  between  these  two 
locks,  and  a  significant  portion  of  the  limestone  foundation.  The  2-D  section  was 
terminated  at  the  center  line  of  the  existing  McAlpine  Lock  and  the  center  line  of  the 
new  lock. 

Backfill  Placement  Method  of  Analysis 

One  of  the  earliest  successful  applications  of  soil-structure  interaction  analysis 
was  performed  by  Clough  and  Duncan  (1969)  in  their  analysis  of  the  two 
reinforced-concrete  U-frame  locks  at  Port  Allen  and  Old  River.  These  two  locks 
had  been  extensively  instrumented.  Prior  to  Clough  and  Duncan's  analysis,  the 
instrumentation  data  had  been  thought  to  be  unreliable  and  contrary  to  the  perceived 
understanding  of  the  behavior  of  locks  to  loadings  encountered  during  lock 
operation.  Clough  and  Duncan's  study  showed  that  the  best  agreement  between 
results  computed  using  the  finite  element  method  and  those  obtained  through 
instrumentation  measurements  is  obtained  when  the  actual  construction  process  is 
simulated  as  closely  as  possible  in  the  analysis.  During  their  study,  Clough  and 
Duncan  developed  what  is  referred  to  as  a  backfill  placement  analysis  in  which  the 
loads  exerted  by  the  backfill  on  the  lock  wall  are  generated  automatically  during 
simulated  placement  of  backfill  behind  the  wall  (i.e.,  predetermined  earth  pressure 
force  distributions  between  the  soil  and  the  lock  are  not  specified).  This  requires 
that  the  soil  backfill  and  foundation  soil  strata  be  included  in  the  finite  element 
mesh.  This  procedure  involved  the  use  of  incremental  finite  element  analysis  with 
nonlinear,  stress-dependent,  stress-strain  behavior  for  the  soil.  Linear  elastic 
behavior  was  assumed  for  the  concrete  lock  wall.  An  additional  requirement  is  that 
interface  elements  be  incorporated  within  the  finite  element  mesh  to  allow  for 
relative  movement  between  the  soil  and  structure.  Since  the  Clough  and  Duncan 
study,  soil-structure  interaction  analysis  using  the  backfill  placement  procedure  has 
been  successfully  applied  to  a  variety  of  earth  retaining  structures.  An  extensive 
database  of  strength  and  hyperbolic  stress-strain  parameters  (using  data  obtained 
from  consolidated-drained  triaxial  tests  on  a  variety  of  soils)  has  also  evolved 
because  of  the  extensive  number  of  applications  of  this  family  of  numerical 
procedures  in  engineering  practice  (see  discussion  in  Ebeling  and  Mosher  1996). 

An  updated  version  of  SOILSTRUCT,  referred  to  as  SOILSTRUCT-ALPHA,  is 
used  in  this  study. 


1  All  elevations  (el)  cited  herein  are  in  feet  referenced  to  lhe  National  Geodetic  Vertical  Datum. 
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Figure  1.  Typical  section,  north  wall  of  new  RCC  McAlpine  Lock  (0.305  m  =  1  ft) 
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Figure  2.  North  wall  of  new  RCC  McAlpine  Lock  (0.305  m  =  1  ft) 


Analysis  Description 

SOILSTRUCT- ALPHA  (Ebeling,  Duncan,  and  Clough  1990)  is  a  special- 
purpose,  finite  element  program  for  2-D,  plane  strain  analysis  of  soil -structure 
interaction  problems.  SOILSTRUCT  calculates  displacements  and  stresses 
resulting  from  incremental  construction,  backfilling,  excavation,  dewatering,  rising 
water  table,  and/or  load  application.  Nonlinear,  stress-path-dependent,  stress-strain 
behavior  of  the  backfill  was  approximated  in  the  finite  element  analysis  using  the 
tangent  modulus  method.  In  the  tangent  modulus  method,  new  values  of  tangent 
moduli  are  assigned  to  each  soil  element  at  each  increment  of  loading  (i.e., 
dewatering,  lock  construction,  and  backfilling)  or  unloading  (i.e.,  excavation,  rising 
water  table).  The  modulus  values  assigned  to  each  element  are  adjusted  in 
accordance  with  their  stresses  to  simulate  nonlinear  behavior. 

SOILSTRUCT  was  expanded  during  the  U.S.  Corps  of  Engineers’  first  Repair, 
Evaluation,  Maintenance,  and  Rehabilitation  (REMR)  Research  Program  to  model 
the  loss  of  contact  between  the  base  of  a  wall  (a  lock  in  this  case)  and  its  rock 
foundation  using  a  procedure  called  the  ALPHA  method  (Ebeling,  Duncan,  and 
Clough  1990;  Ebeling  et  al.  1992).  The  ALPHA  method  was  extended  to  soil 
elements  by  Regalado,  Duncan,  and  Clough  (1992)  to  reduce  numerical  inaccuracies 
in  soil  elements  that  are  at  or  near  failure. 

Another  enhancement  contained  within  this  version  of  SOILSTRUCT  is  the 
reintroduction  of  the  hyperbolic  shear-stress  displacement  relationship  for  the 
interface  element.  Although  present  in  the  original  version  of  SOILSTRUCT 
(Clough  and  Duncan  1969),  this  nonlinear  relationship  was  missing  from  the 
version  of  computer  code  in  which  the  alpha  method  was  incorporated. 

The  continua  elements  used  to  model  the  soil  and  the  soil-to-slructure  interface 
elements  which  may  have  failed  in  shear  at  one  stage  of  loading  have  the  ability  to 
recover  their  shear  stiffness  and  shearing  resistance  as  a  result  of  an  increase  in 
confining  pressures  at  some  later  stage  of  loading  in  this  version  of  SOILSTRUCT- 
ALPHA.  Several  other  improvements  have  been  made  to  the  material  models  and  to 
the  numerical  procedures  implemented  within  SOILSTRUCT- ALPHA  based  on 
experience  gained  at  the  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES) 
in  conducting  SSI  analyses  of  different  types  of  Corps  structures. 

In  summary,  the  ALPHA  version  of  SOILSTRUCT  used  on  this  project  contains 
numerous  improvements  to  the  analytical  procedures  used  for  modeling  various 
aspects  of  features  impacting  the  SSI  of  lock  walls.  Several  of  these  features  are 
deemed  by  the  authors  of  this  report  to  be  critical  to  an  accurate  assessment  of  the 
SSI  of  the  north  wall  section  of  the  new  McAlpine  Lock. 
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Soil-Structure-Foundation  Interaction 
Analysis  of  New  McAipine  Lock 
North  Wall  Section 

Lock  construction,  backfilling,  and  submergence  model 

Project  engineers  at  the  Louisville  District  established  that  the  probable 
scheduling  of  construction  of  the  new  lock  is  placement  of  a  1.5-ft-  (0.46-m-) 
thick  layer  of  RCC  followed  by  placement  of  backfill  to  the  top  of  the  newly 
placed  RCC  lift.  Each  layer  of  RCC  possesses  sufficient  stiffness  at  the  time  of 
placement  to  support  construction  equipment.  This  staged  (incremental) 
construction  proceeds  until  the  final  elevations  (Figure  3)  are  reached  on  the  new 
lock  and  backfill.  Incremental  construction  of  the  RCC  and  placement  of  the 
backfill  behind  the  wall  is  staggered  by  one  lift  in  the  SOILSTRUCT-ALPHA 
incremental  construction/backfill  placement  analysis.  This  altered  construction 
schedule  in  the  analysis  results  in  a  “stiffer”  layer  of  RCC  elements  being  in 
place,  adjacent  to  the  newly  placed  layer  of  soil.  Thus,  the  conditions  in  the  field 
are  modeled  in  the  finite  element  analysis.  The  stages  of  loading  used  in  the 
SOILSTRUCT-ALPHA  analysis  are  described  in  the  following  paragraphs. 

The  complete  SSI  of  the  north  wall  of  the  new  McAipine  Lock  cross  section 
(Figure  2)  is  modeled  in  three  phases.  The  first  phase  of  the  analysis  introduces 
the  self-weight  of  the  existing  locks  and  existing  backfill  to  the  shale  and 
limestone  foundation.  This  is  accomplished  using  the  gravity  turn-on  option  in 
SOILSTRUCT-ALPHA  prior  to  construction  of  the  new  lock  and  placement  of 
new  backfill.  This  initial  phase  is  illustrated  in  Figure  4a  for  load  case  “I.”  The 
second  phase  models  the  incremental  construction  of  the  RCC  lock  wall  and  the 
incremental  placement  of  the  (moist)  backfill  behind  the  wall.  The  third  phase  of 
the  SSI  analysis  models  the  postconstruction,  partial  submergence  of  the  site. 
Table  1  lists  the  initial  stage  of  loading  and  each  of  the  subsequent  38  stages  of 
loading  in  the  SSI  analysis  of  the  north  wall  section  of  the  new  McAipine  Lock. 

The  incremental  construction  of  the  RCC  lock  wall  is  modeled  in  the  first 
24  load  cases  (following  the  gravity  turn-on  analysis),  with  the  incremental 
placement  of  the  moist  backfill  during  22  of  these  load  increments.  Figure  3 
shows  the  thickness  of  the  layers  of  RCC  and  backfill  that  are  used  in  this  second 
phase  of  the  analysis.  The  thickness  of  the  layers  range  from  1.5  ft  (0.46  m)  thick 
at  the  bottom  to  5  ft  (1.5  m)  thick  at  the  top.  Placement  of  the  first  1.5-ft- 
(0.46-m)  thick  layer  of  RCC  of  the  new  lock  wall  precedes  the  placement  of  the 
first  layer  of  backfill  in  the  SOILSTRUCT-ALPHA  analysis,  labeled  load  case  1 
in  Table  1  and  in  Figure  4a.  This  allows  the  first  layer  of  RCC  to  “stiffen”  prior 
to  the  placement  of  the  adjacent  backfill  layer  in  the  analysis.  The  second  load 
case  models  placement  of  the  first  1.5-ft-thick  layer  of  moist  backfill,  labeled  load 
case  2  in  Table  1  and  shown  in  Figure  4a.  The  second,  “soft”  1.5-ft-  (0.46-m) 
thick  layer  of  RCC  is  also  placed  to  el  373  in  load  case  2  to  expedite  the  analysis 
(and  without  compromising  accuracy).  Incremental  construction  of  the  next  layer 
of  “softer”  RCC  layer  precedes  the  placement  of  the  next  layer  of  backfill  by  one 
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igure  3.  Incremental  construction  of  RCC  lock  wail  and  incremental  placement  of  backfill  at  McAlpine  Lock  (0.305 


Figure  4.  Select  load  cases  modeled  in  the  finite  element  analyses 
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lift  throughout  the  analysis  (as  shown  in  Figures  4a  and  4b,  and  Table  1).  This 
allows  SOILSTRUCT- ALPHA  to  account  for  the  stiffer  RCC  immediately 
adjacent  to  the  newly  placed  layer  of  backfill.  Incremental  placement  of  the 
backfill  to  el  425  concludes  with  load  case  23,  and  incremental  construction  of 
the  lock  concludes  with  load  case  24  (Figure  4b). 

The  third  phase  of  SSI  analysis  models  the  postconstruction,  partial  submerg¬ 
ence  of  the  site  and  is  accomplished  using  SOILSTRUCT- ALPHA  in  14  load 
cases.  The  water  table  in  the  existing  and  new  lock  chambers  is  raised 
incrementally  to  el  383,  and  the  water  table  in  the  new  backfill  is  raised 
incrementally  to  el  395  (listed  as  load  cases  25  through  38  in  Table  1).  A 
hydrostatic  water  table  is  assumed  for  the  backfill  in  all  load  cases.  Loadings 
consist  of  buoyancy  forces  acting  within  the  newly  placed  backfill  and  boundary 
water  pressures  acting  normal  to  both  faces  and  base  of  the  new  RCC  lock,  the 
limestone  foundation,  and  the  faces  of  the  two  existing  lock  walls.  A  linear 
(uplift)  pore  water  pressure  distribution  is  applied  along  the  new  RCC  lock-to- 
limestone  foundation  interface,  with  the  uplift  (pressure)  heads  set  equal  to  13  ft 
(3.96  m)  below  the  toe  and  25  ft  (7.6  m)  below  the  heel  of  the  wall. 

Figure  4c  shows  the  incremental  submergence  of  the  backfill  for  load  cases  25 
and  26  and  for  the  final  load  case  of  38.  The  partial  submergence  of  the  site  is 
accomplished  in  the  finite  element  model  using  two  features  in  SOILSTRUCT- 
ALPHA:  application  of  buoyancy  forces  to  the  submerged  backfill  and 
application  of  boundary  water  pressures  along  the  faces  of  impermeable  materials 
(i.e.,  all  locks  and  the  foundation).  The  buoyancy  loads  acting  on  the  submerged 
soil  elements  are  applied  using  the  SEEP  option  in  SOILSTRUCT.  This 
subroutine  is  modified  for  this  analysis  to  account  for  the  difference  between  the 
submerged  and  moist  unit  weights  for  the  soil.  The  water  table  is  raised  in  the 
newly  placed  backfill  in  the  last  14  load  increments  (Table  1).  Water  pressures 
acting  normal  to  both  faces  and  base  of  the  new  RCC  lock,  the  limestone 
foundation,  and  the  faces  of  the  two  existing  lock  walls  are  applied  incrementally 
in  load  cases  32,  35,  and  38,  as  shown  in  Figure  5.  The  interface  elements  used 
to  model  the  RCC  lock-to-limestone  interface  allow  for  the  application  of  pore 
water  pressures  directed  upward  along  the  base  of  the  new  RCC  lock  and 
boundary  water  pressures  directed  downward  along  the  top  face  of  the  limestone 
interface  region.  This  distribution  of  pore  water  pressure  corresponds  to  the 
nonsite-specific  uplift  pressure  distributions  recommended  for  the  design  of  new 
locks  (Headquarters,  Department  of  the  Army  1995)  and  is  consistent  with 
steady-state  confined  flow  along  a  joint  of  constant  aperture  (Ebeling  and  Pace 
1996,  1997). 


Finite  element  mesh 

Figure  6  shows  the  finite  element  mesh  of  the  Figure  2  cross  section  of  the 
north  wall  of  the  new  lock.  The  mesh  comprises  3,784  nodal  points  and 
3,634  continua  and  162  interface  elements.  Of  the  3,634  continua  elements, 
468  model  new  RCC  lock  and  concrete  facing,  2,154  model  the  limestone  and 
shale  foundation,  and  550  model  the  new  backfill.  The  remaining  462  continua 
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Figure  5.  Incremental  boundary  water  pressures  applied  in  three  load  cases 
(0.305  m  =  1  ft) 
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Figure  6.  Finite  element  mesh  of  north  wall  section  of  new  RCC  McAlpine  Lock  (0.305  m  =  1  ft) 


elements  represent  the  two  existing  locks  and  the  existing  backfill  between  these 
locks.  Of  the  162  interface  elements,  30  model  the  new  RCC  lock-to-limestone 
foundation  interface,  25  model  the  limestone-to-new  backfill  interface,  22  model 
the  new  RCC  lock-to-new  backfill  interface,  and  22  model  the  existing  lock-to- 
new  backfill  interface.  Forty-six  of  the  remaining  interface  elements  are  in  the 
limestone  foundation  along  two  vertical  planes,  each  extending  through  the 
limestone  foundation  from  the  heel  of  the  new  RCC  lock  wall  and  from  the 
existing  lock  wall,  respectively.  These  46  interface  elements  are  “locked” 
together  with  the  assignment  of  high  normal  and  shear  stiffnesses.  The  remaining 
17  interface  elements  extend  along  el  370  from  the  front  face  of  the  existing 
center  lock  wall,  to  the  culvert  of  the  existing,  operational  lock  wall  and  are 
locked  together  with  high  normal  and  shear  stiffnesses. 


Material  properties 

Table  2  lists  the  values  assigned  to  the  elastic  parameters  of  the  RCC,  the 
concrete  facing  for  the  new  lock,  the  concrete  of  the  existing  locks,  and  the 
limestone  and  shale  foundation.  The  elastic  parameters  assigned  to  the  RCC  are 
based  upon  data  obtained  for  the  RCC  mixture  used  on  the  Willow  Creek  Dam 
project  and  provided  by  Portland  District  personnel.  The  soil  that  comprises  the 
backfill  for  the  new  RCC  lock  will  be  a  well-compacted,  well-graded,  dense 
granular  backfill.  Moist  and  saturated  unit  weights  equal  126  pcf 
(2,018.33  kg/m3)  and  130  pcf  (2,082.39  kg/m3),  respectively.  Site-specific 
triaxial  test  data  are  unavailable  for  the  backfill.  Material  parameters  are  assigned 
in  the  finite  element  analysis  based  on  empirical  correlations  to  the  results  for 
similar  types  of  soils  (and  with  the  same  density)  for  which  hyperbolic  stress- 
strain  curve  material  parameters  are  available  (e.g.,  Duncan  et  al.  1978).  An 
additional  requirement  for  the  soil  model  was  that  the  assigned  soil  properties 
correspond  to  an  at-rest  earth  pressure  coefficient  K<,  equal  to  0.44,  based  on  the 
Jaky  equation  for  K0  [=1  -  sin  <{>']  with  4>'  equal  to  34  deg.  Appendix  A  describes 
the  calculations  made  that  result  in  the  assignment  of  values  for  the  hyperbolic 
stress-strain  soil  model  of  the  backfill  given  in  Table  3.  These  calculations 
include  a  settlement  analysis  made  using  SOILSTRUCT- ALPHA  of  a  55-ft- 
(16.76-m)  high,  one-dimensional  (1-D)  soil  column  due  to  self-weight  of  the 
moist  soil. 

Similarly,  no  specific  tests  are  performed  to  define  the  hyperbolic  shear  stress- 
relative  displacement  relationship  for  the  interface  element  used  in 
SOILSTRUCT-ALPHA.  The  assignment  of  material  parameters  to  the  interface 
elements  are  based  on  empirical  correlations  to  interfaces  having  similar  types  of 
soils  (and  with  the  same  density)  for  which  hyperbolic  stress-strain  curve  material 
parameters  are  available  (e.g.,  Clough  and  Duncan  1969,  and  Peterson  et  al. 

1976).  Table  4  summarizes  the  interface  model  and  strength  parameters  assigned 
to  the  RCC  lock-to-limestone  foundation.  Zero  tensile  strength  is  assumed  for  the 
material  comprising  the  RCC  lock-to-limestone  foundation  interface  in  this 
analysis.  Table  5  summarizes  the  material  properties  assigned  to  the  interfaces 
between  the  RCC-to-backfill,  the  limestone-to-backfill,  and  the  existing  lock  wall 
concrete-to-backfill  regions. 
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Table  2 

Elastic  Material  Properties  for  2-D  Elements  Comprising 
SOILSTRUCT-ALPHA  Finite  Element  Model  of  New  RCC  McAlpine 
Lock  (0.517  kN/m3  =  1  pcf,  6.894  psi  =  1  kPa) 

Material  Type 

Unit  Weight,  pcf 

E,  psi 

V 

RCC 

150 

2,000,000 

0.21 

Concrete  Face 

145 

3,400,000 

0.15 

Existing  Concrete 

145 

3,400,000 

0.15 

Limestone 

- 

5,000,000 

0.25 

Shale 

-- 

3,000,000 

0.25 

During  the  postconstruction  partial  submergence  of  the  backfill  to  el  395 
(identified  as  load  cases  25  through  38  in  Table  1),  the  buoyancy  forces  act 
upward,  unloading  the  backfill.  This  results  in  a  reversal  in  the  direction  of  the 
applied  shear  stress  increment  or,  equivalently,  unloading  of  both  the  new  RCC 
lock-to-backfill  interface  elements  and  the  existing  lock  concrete-to-backfill 
interface  elements.  Interface  tests  reported  in  Peterson  et  al.  (1976)  show  that  the 
interface  elements  follow  a  steeper  shear  stress  versus  relative  shear  displacement 
curve  than  that  described  by  the  “primary”  hyperbolic  loading  curve  (the  tangent 
shear  stiffness  kst  =  ksi[I  -  RfiSL]2 ).  To  model  this  behavior,  the  interface  model 
parameter  Rfi  is  set  equal  to  zero  for  these  interface  elements  during  unloading. 
This  eliminates  the  “shear  softening”  while  maintaining  the  influence  of  the 
effective  normal  stress  on  the  value  of  the  tangent  shear  stiffness  k*  by  means  of 
the  term  ksi  (Table  4).  This  concept  is  discussed  further  in  Appendix  B. 


Results  of  Soil-Structure-Foundation 
Interaction  Analysis  of  New  McAlpine 
Lock  North  Wall  Section 

Following  the  gravity  tum-on  analysis,  the  incremental  construction  of  the 
RCC  lock  wall  and  placement  of  moist  backfill  is  modeled  in  the  first  24  load 
cases  for  the  finite  element  model  of  the  McAlpine  lock  shown  in  Figure  6.  This 
is  followed  by  the  postconstruction,  partial  submergence  of  the  site  in  14  load 
cases.  The  results  of  this  complete  SSI  analysis  follow. 


Effective  normal  and  shear  stress  distributions 
computed  at  base  of  new  lock  wall 

Figures  7  and  8  show  the  resulting  effective  normal  and  shear  stress  distribu¬ 
tions  along  the  RCC  lock-to-limestone  interface  after  construction  of  the  new  lock 
and  backfilling  to  el  425  (load  case  24).  Both  effective  normal  and  shear 
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Table  3 


Hyperbolic  Stress-Strain  and  Strength  Parameters  for  Engineered  Backfill 
(0.157  kN/m3  =  1  pcf,  47.88  Pa  =  1  psf) 


Strength 

Parameters 

Hyperbolic  Parameters 

Backfill 

Unit  Weight, 
psf 

c’,  psf 

<t>\  deg 

K. 

K 

n 

kur 

Ks 

m 

Moist 

Granular 

Backfill 

126 

0 

34 

0.44 

500 

1000 

175 

0.025 

Submerged 

Granular 

Backfill 

130 

0 

34 

0.44 

500 

1000 

175 

0.025 

Note: 


Kc  by  Jaky  =  1  -  sin  4>' 

Tangent  Modulus,  E,  =  £,  (  1  -  RF  SL  )2 

'  V 

Initial  Modulus,  E,  =  K  P„  — 

P‘  , 

Stress  Level,  SL  «  (  o,  -  o3  )  /  (  a,  -  o3  )Mur. 


(  °1  *  °3  )  Failure  “ 


2  c  cos  4>  +  2  Og  sin  4> 
1  -  sin  4> 


Unload  Reload  Modulus,  EUR  =  KUR  Pa 


V  P‘> 


Bulk  Modulus,  B  = 


(3  -  6  v_) 


Bulk  Modulus,  B  =  KB  P,  ^  —  j 

vnom  =  Nominal  value  of  Poissorfs  ratio 

Poissorfs  ratio,  v  =  -1  |  1  -  [  (  1  -  2  )  (  1  -  RF  SL  f  ]  ] 

P*  -  atmospheric  pressure 

computed  below  the  toe  of  the  new  RCC  lock  (x  =  0  ft)  because  the  greatest  con¬ 
crete  mass  is  concentrated  above  this  region  and  because  of  the  overturning  moment 
about  the  toe  resulting  from  the  55  ft  (16.76  m)  of  backfill.  The  effective  normal 
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Table  4 

Material  Properties  for  the  RCC-to-Limestone  Interface  Elements 
of  the  SOILSTRUCT-ALPHA  Finite  Element  Model  of  McAlpine 
Lock  (47.88  Pa  =  1  psf,  1  MPa/m  =  3.684  psi/in.,  1  MPa/m  = 
6,365.9  psf/ft) 


Strength  Parameters 

Stiffness  Parameters 

Interface  Region 

c',  psf 

deg 

k,,  psi/in. 

K,  psf/ft 

RCC-to-Limestone 

2300 

41 

10,000 

1 .0  x  1 0s 

Equations  for  Interface  Model 


The  normal  stress  at  the  center  of  the  interface  element  is  given  by 

=  kn 

where  A„  is  the  average  relative  displacement  normal  to  the  interface  element.  The  shear  stress  at 
the  center  of  the  interface  element  is  given  by 

T,  ■  K  A, 

where  At  is  the  average  relative  shear  displacement  along  the  interface  element. 


effective  normal  stress  below  the  heel  of  the  new  RCC  lock  (x  =  47  ft)  is  equal  to 
1,626  psf  (7,939  kg/m2).  Flooding  of  the  lock  chamber  to  el  383  and  partial 
submergence  of  the  backfill  to  el  395  (load  case  38)  reduce  the  effective  normal 
stresses  and  increase  the  shear  stresses  along  the  base,  as  shown  in  Figures  7  and 
8.  The  effective  normal  stress  below  the  heel  of  the  new  RCC  lock  is  reduced  to 
0  psf  (kg/m2);  however,  full  base  area  contact  is  maintained. 

Results  after  backfilling  to  el  425 

Figures  9  and  10  show  the  variation  of  horizontal  effective  stress  o\  and  shear 
stress  xxy  with  elevation  along  a  vertical  plane  extending  through  the  backfill  from 
the  heel  of  the  new  RCC  lock  wall  (x  =  47  ft  (14.33  m))  after  backfilling  to  el  425 
(load  case  24  in  Table  1).  This  section  is  designated  as  section  A-A  (see 
Figure  11)  throughout  this  report.  As  anticipated,  these  figures  show  the  resulting 
effective  horizontal  and  shear  stress  distributions  to  increase  with  depth  below 
the  surface  of  the  backfill. 

The  resultant  horizontal  effective  force  of  the  Figure  9  distribution  of  o\ ,  F„ , 
is  equal  to  83,367  lb  per  ft  run  of  wall  (1,216.7  kN  per  m)  and  acts  at  el  390.28 
(=  0.37HtackM)  along  section  A-A.  The  horizontal  earth  pressure  coefficient  Kj,  is 
computed  along  section  A-A  using  the  relationship 


Kh 


Effective  Overburden 


(1) 
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Table  5  .  . 

Material  Properties  for  Interface  Elements  Comprising 
SOILSTRUCT-ALPHA  Finite  Element  Model  of  RCC  McAlipine 
Lock  for  Backfill  Placement  Analysis  (1  MPa/m  =  6,365.9  psffft) 


Hyperbolic  Parameters 


Interface  Region 

♦'i.  deg 

n. 

Rs 

RCC-to-Backfill 

34 

1.0x10“ 

0.8 

0.9 

Concrete-to-Backfill 

34 

1.0  x  10“ 

0.8 

0.9 

Limestone-to-Backfill 

34 

1.0  x  10“ 

0.8 

0.9 

Normal 

Stiffness 

k„,  psf/ft 

1.0  x  108 
1.0  x  108 
1.0  x  10® 


Equations  for  the  Interface  Model 

The  normal  stress  at  the  center  of  the  interface  element  is  given  by 


=  kn  A„ 


where  A„  is  the  average  relative  displacement  normal  to  the  interface  element.  For  each  load 
increment,  the  change  in  shear  stress  at  the  center  of  the  interface  element  is  given  by 


&ts  =  ksl  As 


where  A,  is  the  average  change  relative  shear  displacement  along  the  interface  element. 
^J!  =  ksi  (  1  Rfi  SL-i ) 


■  Kl  pT  j 

SL,  -  T  /  TFa(/ura  =  T  /  (  o'  tan  <t)( ) 
yw  =  unit  weight  of  water  P 


PA  =  atmospheric  pressure 


where  the  Effective  Overburden  is  defined  as 

rEl  425  /  » 

Effective  Overburden  =  j  ooverburden  ay  (2) 

J  El  370 

In  the  case  of  a  55-ft-  (16.76-m)  high  column  of  moist  backfill  with  equal  to 
126  pcf  (2,018.3  kg/m3),  Equation  2  for  the  Effective  Overburden  becomes 


Effective  Overburden  =  —  (H) 
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Figure  7.  Effective  normal  stress  distributions  along  the  base  of  new  RCC  lock  wall  (0.305  m  =  1  ft, 
47.88  Pa  =  1  psf) 

and  is  equal  to  190,575  lb  per  ft  run  of  wall  (2,781.2  kN  per  m).  The  value  for  K„ 
after  backfilling  to  el  425  is  equal  to  0.437  by  Equation  1. 

Calculation  of  the  resultant  horizontal  effective  force  Fx  of  the  resulting 
distributions  of  o'x  with  elevation  in  the  backfill  is  made  after  backfilling  to  el  425 
for  two  other  vertical  sections.  Figure  1 1  shows  the  locations  of  the  three  sections 
used  in  the  summary  of  the  results.  The  two  additional  sections  are  designated  as 
section  B-B  at  x  =  74.5  ft  (22.7  m),  and  section  C-C  at  x  =  100.9  ft  (30.75  m). 

Table  6  summarizes  the  results  of  calculations  made  to  determine  the 


22 


Soil-Structure-Foundation  Interaction  Analysis 


o 

<D 

_C 

CO 


Figure  8.  Shear  stress  distributions  along  the  base  of  new  RCC  lock  wall  (0.305  m  =  1  ft,  47.88  Pa  = 
1  psf) 

resultant  horizontal  effective  force  Fx  and  horizontal  earth  pressure  coefficient  K,,  at 
sections  A-A,  B-B,  and  C-C.  The  table  shows  the  resultant  horizontal  effective 
force  Fx  for  these  two  distributions  of  a'x  (not  shown)  result  in  values  of  F„  equal  to 
83,520  lb  per  ft  run  of  wall  (1,218.9  kN  per  m)  (section  B-B)  and  71,198  lb  per  ft 
run  of  wall  (1,039.1  kN  per  m)  (section  C-C).  Figure  1 1  also  shows  the  variation  in 
horizontal  earth  pressure  coefficient  K,, ,  within  distance  from  the  back  of  the  new 
RCC  lock  wall  for  the  three  sections.  Figure  1 1  shows  the  values  for  Kj,  equals 
0.437  adjacent  to  the  new  RCC  lock  wall  (section  A-A,  0.438  at  section  B-B,  and 
0.374  adjacent  to  the  existing  lock  wall  (section  C-C). 
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Figure  9.  Horizontal  effective  stress  distributions  along  section  A-A  after  backfilling  to  el  425  (0.305  m  = 
1  ft,  47.88  Pa  =  1  psf) 


Note  that  the  values  for  at  sections  A-A  and  B-B  are  just  slightly  less  than  0.44, 
the  value  for  K0  from  the  1-D  soil  column  settlement  analysis  reported  in  Appendix 
A.  Recall  that  the  conditions  corresponding  to  a  Ko  stress  state  within  the  soil  exist 
in  a  region  of  constrained,  uniform  settlement.  The  value  of  K*,  at  section  C-C  is  1 8 
percent  less  than  the  K0  value  because  of  the  influence  of  the  large  vertical  shear 
force  acting  along  this  plane  (to  be  discussed  subsequently)  resulting  from  the 
presence  of  the  existing  concrete  lock  wall  at  x  =  102  ft  (31.09  m). 

The  resultant  vertical  shear  force  of  the  Figure  10  distribution  of  ,  Fv ,  is  equal 
to  1 1,416  lb  per  ft  run  of  wall  (166.6  kN  per  m)  along  section  A-A.  The 
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Elevation  (ft) 


Figure  10.  Shear  stress  distributions  along  section  A-A  after  backfilling  to  el  425  (0.305  m  =  1  ft, 
47.88  Pa  =  psf) 


vertical  earth  pressure  coefficient  K,  is  computed  along  section  A-A  using  the 
relationship 


Effective  Overburden  W 

where  the  Effective  Overburden  is  equal  to  190,575  lb  per  ft  run  of  wall 
(2,781.2  kN  per  m)  by  Equation  3.  The  value  for  K,  after  backfilling  to  el  425  is 
equal  to  0.06  by  Equation  4. 
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Elevation  (ft) 


Figure  1 1 .  Variation  of  horizontal  earth  pressure  coefficient,  K„  with  distance  from  heel  of  new  RCC 
lock  wall  after  backfilling  to  el  425  (0.305  m  =  1  ft) 
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Table  6 

Summary  of  SOILSTRUCT-ALPHA  Results  After  Backfilling  New 
RCC  McAlpine  Lock  to  el  425  (14.594  N  per  m  =  1  lb  per  ft  run  of 
wall) 

Section 

Effective  Overburden, 
lb  per  ft  run  of  wall 

F,‘,  lb  per  ft 
run  of  wall 

Fv,  lb  per  ft 
run  of  wall 

1C 

K, 

A-A 

190,575 

83,367 

11,416 

0.437 

0.060 

B-B 

190,575 

83,520 

-4,990 

0.438 

-0.026 

C-C 

190,575 

71,198 

-34,869 

0.374 

-0.183 

Calculation  of  the  resultant  vertical  shear  force  Fv  from  the  resulting  distribu¬ 
tions  of  Txy  with  elevation  is  made  after  backfilling  to  el  425  for  sections  B-B  and 
C-C.  Table  6  summarizes  the  results  of  calculations  made  to  determine  the 
resultant  vertical  shear  force  Fv  and  vertical  earth  pressure  coefficient  K*  at 
sections  A-A,  B-B,  and  C-C.  The  table  shows  that  the  resultant  vertical  shear 
force,  Fv ,  for  these  two  distributions  of  txy  (not  shown)  result  in  values  of  Fv 
equal  to  4,990  lb  per  ft  run  of  wall  (72.8  kN  per  m)  (section  B-B)  and  34,869  lb 
per  ft  run  of  wall  (508.9  kN  per  m)  (section  C-C).  Figure  12  shows  the  variation 
in  vertical  earth  pressure  coefficient  1C,  within  distance  from  the  back  of  the  new 
RCC  lock  wall  for  the  three  sections.  Figure  12  shows  the  values  for  K,,  equal 
0.06  adjacent  to  the  new  RCC  lock  wall  (section  A-A),  0.026  at  section  B-B,  and 
0.183  adjacent  to  the  existing  lock  wall  (section  C-C).  The  values  for  Ky  for 
sections  B-B  and  C-C  are  plotted  as  negative  values  to  note  that  the 
corresponding  vertical  shear  forces  at  these  sections  act  counter  to  the  direction  of 
Fv  along  section  A-A,  as  shown  in  the  three  diagrams  in  Figure  12.  Note  that  had 
there  been  no  differential  settlement  across  a  vertical  section  in  the  backfill,  Fv 
and  Ky  would  be  equal  to  0.  (This  is  the  case  for  the  settlement  of  the  Appendix 
A  1-D  soil  column  due  to  self-weight.)  Fv  and  Ky  are  nonzero  at  the  three  sections 
because  of  the  close  proximity  of  the  two  shear  faces  furnished  by  the  new  RCC 
lock  and  the  existing  lock  walls. 

Table  7  summarizes  the  resultant  forces  and  their  points  of  application  on  the 
free  body  of  the  new  RCC  lock  wall  shown  in  Figure  13  after  backfilling  to 
el  425.  These  resultant  forces  and  their  points  of  application  are  computed  by 
SOILSTRUCT-ALPHA  using  the  distributions  of  effective  normal  and  shear 
stresses  computed  within  the  interface  elements  located  along  the  back  of  the  new 
RCC  lock  wall  and  along  the  base  of  the  lock  wall  for  load  case  24.  Twenty-two 
interface  elements  delineate  the  back  of  the  RCC  lock  wall,  and  thirty  interface 
elements  delineate  the  base  of  the  lock  wall.  The  Table  7  values  of  forces  and 
their  points  of  application  are  consistent  with  the  actual  distribution  of  stresses 
along  the  two  planes  defined  by  the  two  groups  of  interface  elements  because  of 
the  analytical  formulation  used  to  define  the  interface  elements  in  SOILSTRUCT- 
ALPHA.  The  Table  7  results  show  that  the  amount  of  shear  (or  “downdrag”) 
along  the  back  of  the  new  RCC  lock  wall  is  significant  and  corresponding  to  a 
mobilized  angle  of  interface  friction  8^  equal  to  25.28  deg. 
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Figure  12.  Variation  of  vertical  earth  pressure  coefficient,  K*  with  distance  from  heel  of  new  RCC  lock 
wall  after  backfilling  to  el  425  (0.305  m  =  1  ft) 
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Table  7 

Summary  of  SOU 
El  425  (14.594  N 

.STRUCT-ALPHA  Results  After  Backfilling  New  RCC  McAlpine  Lock  to 
per  m  =  1  lb  per  ft  run  of  wall,  0.305  m  =  1  ft) 

Load 

Case 

F„  lb  per  ft 
run  of  wall 

Fn,  lb  per  ft  run 
of  wall 

WL 

WL 

deg 

T,  lb  per  ft 
run  of  wall 

N',  lb  per  ft 
run  of  wall 

Xjtt 

&..B 

24 

66,435 

140,661 

21.96 

0.34 

25.28 

83,003 

405,573 

17.98 

0.38 

Note: 

where  B  =  47  ft,  L  =  65.19  ft,  and  tan  (  bmob  )  =  — 

F' 

Results  after  partial  submergence  to  el  395 

Figures  14  and  15  show  the  variation  of  horizontal  effective  stress  o',  and  shear 
stress  Txy  with  elevation  along  section  A-A  after  backfilling  to  el  425  and  after 
partial  submergence  of  the  backfill  to  el  395  (load  cases  24  and  38,  respectively). 

As  anticipated,  these  figures  show  that  the  resulting  effective  horizontal  and  shear 
stress  distributions  increase  with  depth  below  the  surface  of  the  backfill.  Figure  14 
shows  that  the  magnitude  of  the  values  of  o\  below  el  395  decreases  with  the  partial 
submergence  of  the  site.  Figure  15  shows  that  the  magnitude  of  the  values  of  Txy 
decreases  with  the  partial  submergence  of  the  site^  During  postconstruction,  after 
partial  submergence  of  the  backfill  to  el  395,  the  buoyancy  forces  act  upward, 
thereby  unloading  the  backfill.  This  also  results  in  a  reversal  in  shear  stress 
increment. 


The  resultant  horizontal  effective  force  of  the  Figure  14  distribution  of  o'x,  Fx,  is 
equal  to  79,102  lb  per  ft  run  of  wall  (1,154.4  kN  per  m)  and  acts  at  el  390.99 
(=  0.381^^1,)  along  section  A-A  after  partial  submergence  of  the  backfill  to 
el  395.  The  horizontal  earth  pressure  coefficient  is  computed  along  section  A-A 
using  Equation  1,  with  the  Effective  Overburden  computed  using  Equation  2.  In  the 
case  of  partially  submerged,  55-ft-  (16.76-m-)  high  column  [with  ymoist  equal  to 
126  pcf  (2,018.32  kg/m3)  and  equal  to  130  pcf  (2,082.39  kg/m3)]. 

Equation  2  for  the  Effective  Overburden  becomes 


=  \  y ( Di  )2  +  ymoin  ( Di  d2  >  +  \  ykuy  ( Di  )2 

where 


(5) 


Dj  =  the  thickness  of  backfill  above  the  hydrostatic  water  table  =  30  ft  (9  m) 

D2  =  the  thickness  of  submerged  backfill  above  the  base  of  the  wall  =  25  ft 
(7.6  m) 
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Figure  13.  Planes  along  surfaces  of  the  new  RCC  lock  on  which  resultant 

forces  are  computed  after  backfill  is  placed  to  el  425  (0.305  m  =  1  ft) 


Ybuy  =  buoyant  unit  weight  of  submerged  backfill,  y  saturated  -  Yw 

yw  =  unit  weight  of  water  =  62.4  pcf  (999.55  kg/m3) 

and  is  equal  to  172,325  lb  per  ft  run  of  wall  (2,5 14.9  kN  per  m).  The  value  for  K„ 
after  partial  submergence  of  the  backfill  to  el  395  is  equal  to  0.459  by  Equation  1. 

Calculation  of  the  resultant  horizontal  effective  force  Fx  of  the  resulting 
distributions  of  o\  with  elevation  in  the  backfill  is  made  after  partial  sub¬ 
mergence  of  the  backfill  to  el  395  for  sections  B-B  and  C-C.  Table  8  summarizes 
the  results  of  calculations  made  to  determine  the  resultant  horizontal  effective 
force  Fx  and  horizontal  earth  pressure  coefficient  K*  at  sections  A-A,  B-B,  and 
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Figure  14.  Horizontal  effective  stress  distributions  along  section  A-A  (0.305  m  =  1  ft,  47.88  Pa  =  1  psf) 


C-C.  The  table  shows  the  resultant  horizontal  effective  force  Fx  for  these  two 
distributions  of  o'x  (not  shown)  result  in  values  of  Fx  equal  to  79,022  lb  per  ft  run 
of  wall  (1,153.2  kN  per  m)  (section  B-B)  and  68,475  lb  per  ft  run  of  wall 
(999.3  kN  per  m)  (section  C-C).  Figure  16  shows  the  variation  in  horizontal  earth 
pressure  coefficient  K*, ,  with  distance  from  the  back  of  the  new  RCC  lock  wall 
for  the  three  sections.  Figure  16  shows  the  values  for  K*  equals  0.459  adjacent  to 
the  new  RCC  lock  wall  (section  A-A),  0.459  at  section  B-B,  and  0.397  adjacent 
to  the  existing  lock  wall  (section  C-C).  The  values  for  K*  after  partial 
submergence  of  the  backfill  to  el  395  (Figure  16)  are  greater  than  the  values  for 
K,,  after  backfilling  to  el  425  (Figure  1 1)  as  a  result  of  the  “unloading”  of  the  soil 
during  submergence. 
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Figure  15.  Shear  stress  distributions  along  section  A-A  (0.305  m  =  1  ft,  47.88  Pa  =  1  psf) 


The  resultant  vertical  shear  force  of  the  Figure  15  distribution  of  x^ ,  Fv ,  is 
equal  to  6,775  lb  per  ft  run  of  wall  (98.9  kN  per  m)  along  section  A-A.  The  vertical 
earth  pressure  coefficient  is  computed  along  section  A-A  using  Equation  4  with 
the  Effective  Overburden  equal  to  172,325  lb  per  ft  run  of  wall  (2,514.9  kN  per  m) 
by  Equation  5.  The  value  for  after  partial  submergence  of  the  backfill  to  el  395 
is  equal  to  0.039  by  Equation  4. 

Calculation  of  the  resultant  vertical  shear  force  Fv  of  the  resulting  distributions 
of  x^  with  elevation  in  the  backfill  is  made  after  partial  submergence  of  the  backfill 
to  el  395  for  sections  B-B  (x  =  74.5  ft  (22.7  m))  and  section  C-C  (x  = 
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Table  8 

Summary  of  SOILSTRUCT-ALPHA  Results  After  Backfilling  New  RCC  McAlpine  Lock 
to  El  425  and  After  Submergence  to  El  395  (14.594  N  per  m  =  1  lb  per  ft  run  of  wall) 

Section 

Effective  Overburden,  lb 
per  ft  run  of  wall 

F,',  lb  per  ft  run  of 
wall 

Fv,  lb  per  ft 
run  of  wall 

K* 

Kv 

A-A 

172,325 

79,102 

6,775 

0.459 

0.039 

B-B 

172,325 

79,022 

-4,241 

0.459 

-0.025 

C-C 

172,325 

68,475 

-25,389 

0.397 

-0.147 

100.9  ft  (30.75  m)),  with  the  results  summarized  in  Table  8.  The  table  shows  the 
resultant  vertical  shear  force  Fv  for  these  two  distributions  of  txy  (not  shown) 
result  in  values  of  Fv  equal  to  4,241  (section  B-B)  and  25,389  (section  C-C). 
Figure  17  shows  the  variation  in  vertical  earth  pressure  coefficient  K„ ,  within 
distance  from  the  back  of  the  new  RCC  lock  wall  for  the  three  sections.  Figure 
17  shows  the  values  for  1C,  equal  0.039  adjacent  to  the  new  RCC  lock  wall 
(section  A- A),  0.025  at  section  B-B,  and  0.147  adjacent  to  the  existing  lock  wall 
(section  C-C).  The  values  for  1C,  after  partial  submergence  of  the  backfill  to  el 
395  (Figure  17)  are  less  than  the  values  for  Kv  after  backfilling  to  el  425 
(Figure  12).  This  behavior  is  a  result  of  the  buoyancy  forces  acting  upward, 
thereby  unloading  the  backfill  and  reversing  the  direction  of  the  applied  shear 
stresses  (as  compared  to  their  direction  during  backfilling)  along  the  faces  of  the 
new  RCC  and  existing  lock  walls. 

Table  9  summarizes  the  resultant  forces  and  their  points  of  application  on  the 
free  body  of  the  new  RCC  lock  wall  through  which  imaginary  section(s)  (shown 
in  Figure  18)  are  made.  These  resultant  forces  and  their  points  of  application  are 
computed  by  SOILSTRUCT-ALPHA  after  backfilling  to  el  425  (load  case  24  in 
Table  1)  and  after  partial  submergence  of  the  backfill  to  el  395  (load  case  38). 

The  Table  9  results  show  that  the  shear  force  Fs  along  the  back  of  the  new  RCC 
lock  wall  is  reduced  by  16.6  percent  with  submergence  of  the  backfill  to  el  395  as 
compared  to  the  results  in  Table  7  (load  case  24).  The  mobilized  angle  of 
interface  friction  6^*  is  reduced  from  25.28  to  23.33  deg  with  partial  submer¬ 
gence  of  the  backfill. 


Summary 

This  report  summarizes  the  results  of  a  complete  soil-structure  interaction 
analysis  of  the  new  RCC  Me  Alpine  Lock  using  SOILSTRUCT-ALPHA.  The 
principal  results  are  as  follows: 

a.  A  complete  soil-structure  interaction  analysis  of  the  new  RCC  McAlpine 
Lock  results  in  nonlinear  effective  normal  and  shear  stress  distributions 
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Figure  16.  Variation  of  horizontal  earth  pressure  coefficient,  K„  with  distance  from  heel  of  new 
RCC  lock  wall  after  submergence  to  el  395  (0.305  m  =  1  ft) 
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Fv  =  Ky  •  (Effective  Overburden) 
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*  Negative  values  of  Kv  reflect 
directed  upward  on  sections 
B  —  B  and  C  —  C 


Figure  1 7.  Variation  of  vertical  earth  pressure  coefficient,  K„  with  distance  from  heel  of  new 
RCC  lock  wall  after  submergence  to  el  395  (0.305  m  =  1  ft) 
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Table  9 

Summary  of  SOILSTRUCT-ALPHA  Results  After  Backfilling  New  RCC  McAlpine  Lock 
to  El  425  and  After  Submergence  to  El  395  (14.594  N  per  m  =  1  lb  per  ft  run  of  wall, 

0.305  m  -  1  ft) 

For  Load  Cases  24  and  38 

Load 

Case 

F.,  lb  per 
ft  run  of 
wall 

F'„,  lb  per 
ft  run  of 
wall 

iFnm 

LA 

deg 

T,  lb  per 
ft  run  of 
wall 

N',  lb  per 
ft  run  of 
wall 

x«,ft 

XJB 

24 

66,435 

21.96 

0.34 

25.28 

83,003 

405,573 

17.98 

0.38 

38 

55,404 

23.32 

0.36 

23.33 

illlll® 

347,902 

16.18 

0.34 

For  Load  Case  38 

Load 

Case 

U^lb 
per  ft 
run  of 
wall 

Xummi  ft 

B 

lb 

per  ft  run 
of  wall 

■ 

U^telb 
per  ft 
run  of 
wall 

lub«*,ft 

38 

25.97 

0.55 

5,272.8 

4.33 

23,113 

9.88 

where  B  =  47  ft,  L  =  65.19  ft,  and  tan  (5^  )  =  — ’ 

along  the  base  of  the  lock  at  the  end  of  backfilling  to  el  425,  and  after 
flooding  the  lock  to  el  383  and  partial  submergence  of  the  backfill  to 
el  395. 

b.  Full  base  contact  is  maintained  throughout  all  (38)  stages  of  loading. 
However,  the  effective  normal  stress  below  the  heel  of  the  new  RCC  lock 
is  reduced  to  0  (kPa)  psf  after  flooding  the  lock  to  el  383  and  partial 
submergence  of  the  backfill  to  el  395. 

c.  Among  the  results  computed  in  a  backfill  placement  analysis  is  the 
distribution  of  horizontal  effective  stress  at  three  vertical  sections  in  the 
backfill.  At  section  A- A,  the  value  of  the  effective  horizontal  force  Fx  is 
computed  from  the  distribution  of  effective  horizontal  stress  o\  along  a 
vertical  plane  extending  through  the  backfill  from  the  heel  of  the  wall. 
The  value  for  Fx  is  also  characterized  in  terms  of  the  value  of  the 
horizontal  earth  pressure  coefficient  K,,.  The  results  of  the  SSI  analysis 
shows  that  Fx  decreases  as  a  result  of  the  submergence  of  the  backfill  and 
the  value  of  Kj,  increases. 

d.  The  value  of  K+,  at  section  C-C  is  less  than  the  values  of  K*  at  sections 
A- A  and  B-B  because  of  the  influence  of  the  large  vertical  shear  force 
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Figure  1 8.  Planes  along  surfaces  of  the  new  RCC  lock  on  which  resultant  forces  are  computed  - 
submergence  to  el  395  (0.305  m  =  1  ft) 


along  this  plane,  caused  by  the  presence  of  the  existing  lock  wall.  Recall 
that  section  C-C  is  immediately  adjacent  to  the  existing  lock  wall. 

e.  The  distribution  of  shear  stress  at  three  vertical  sections  in  the  backfill  is 
also  computed  in  a  backfill  placement  analysis.  At  section  A-A,  the  value  of 
the  vertical  shear  force  Fv  is  computed  from  the  distribution  of  shear  stress 
along  a  vertical  plane  extending  through  the  backfill  from  the  heel  of  the 
wall.  The  value  for  Fv  is  also  characterized  in  terms  of  the  value  of  the 
vertical  earth  pressure  coefficient  The  results  of  the  SSI  analysis  shows 

that  Fv,  or,  equivalently,  K^,  reduces  as  a  result  of  the  submergence  of  the 
backfill. 
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/.  This  SSI  analysis  shows  that  in  the  case  of  a  postconstruction  rise  in  the 
groundwater  level  in  the  backfill,  a  rebound  of  the  soil  occurs.  This  results 
in  a  reduction  in  effective  stress  in  the  backfill  and  a  reduction  in  the  shear 
force  Fv.  This  also  occurred  at  Red  River  Lock  No.  1  (Ebeling  et  al.  1993, 
or  Ebeling  and  Mosher  1996). 

g.  During  the  postconstruction,  partial  submergence  of  the  backfill  to  el  395, 
the  buoyancy  forces  act  upward,  unloading  the  backfill.  This  results  in  a 
reversal  in  the  direction  of  the  applied  shear  stress  increment  or, 
equivalently,  unloading  of  both  the  new  RCC  lock-to-backfill  interface 
elements  and  the  existing  lock  concrete-to-backfill  interface  elements. 
Interface  tests  reported  in  Peterson  et  al.  (1976)  show  that  the  interface 
elements  follow  a  steeper  shear  stress  versus  relative  shear  displacement 
curve  than  that  described  by  the  “primary”  hyperbolic  loading  curve. 
Incorporating  this  type  of  interface  behavior  in  the  SSI  analysis  of  the  new 
RCC  McAlpine  Lock  wall  is  important  to  the  accuracy  of  the  computed 
results. 

h.  Appendix  B  reports  on  the  errors  related  to  using  the  “softer”  primary 
hyperbolic  curve  for  the  interface  elements  during  partial  submergence  of 
the  backfill  to  el  395  compared  to  the  results  from  using  the  more 
appropriate  “stiffer”  shear  stress  versus  relative  shear  displacement  curve. 
The  most  important  error  resulting  from  using  a  soft  shear  stress  versus 
relative  shear  displacement  curve  during  unloading  for  the  interface  elements 
is  in  the  overprediction  of  the  value  for  the  effective  normal  stress  computed 
below  the  heel  of  the  new  RCC  lock  wall. 
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Appendix  A 

Backfill  Placement  Analysis  of 
Moist  One-Dimensional  Soil 
Column 


There  are  two  major  requirements  for  the  assignment  of  material  parameters 
in  the  complete  soil-structure  interaction  analysis  of  the  new  roller-compacted 
concrete  (RCC)  McAlpine  Lock  wall  using  the  backfill  placement  method  that  is 
incorporated  in  SOILSTRUCT-ALPHA.  One  requirement  is  that  the  stress- 
strain  model  for  the  soil  representing  the  backfill  be  representative  of  the  soil 
comprising  the  backfill.  The  second  requirement  is  that  the  material  properties 
assigned  to  the  soil  used  to  model  the  backfill  possess  an  at-rest  earth  pressure 
coefficient  equal  to  0.44  (from  Jaky’s  relationship  for  K0  (=  1  -  sin  <j)'>  with  <J)' 
equal  to  34  deg).  Calculations  described  in  this  appendix  are  used  in  the 
assignment  of  the  value  of  the  nominal  Poisson's  ratio  for  soil. 

This  appendix  summarizes  the  results  of  the  backfill  placement  analysis  of  the 
moist  one-dimensional  (1-D)  soil  column  shown  in  Figure  Al.  Settlement  of  the 
soil  column  resulted  from  self-weight.  Calculations  are  made  using 
SOILSTRUCT-ALPHA. 

Site-specific  biaxial  test  data  were  unavailable  for  the  new  RCC  lock  wall 
backfill.  Material  parameters  are  assigned  in  the  finite  element  analysis  based  on 
empirical  correlations  to  the  results  for  similar  types  of  soils  (and  with  the  same 
density)  for  which  hyperbolic  stress-strain  curve  material  parameters  are  avail¬ 
able  (e.g.,  Duncan  et  al.  1978  ‘).  The  values  listed  in  Table  Al  are  the  result  of 
the  evaluation  made  for  this  project. 

A  free-field  soil  column  is  constructed  to  be  consistent  with  the  layering  used 
in  the  Figure  6  (main  text)  two-dimensional  finite  element  model  for  the  new 
backfill.  The  compression  (settlement)  of  the  Figure  Al  1-D  (constrained)  soil 
column  has  the  following  characteristics: 


1  References  cited  in  this  appendix  are  listed  in  the  References  at  the  end  of  the  main  text. 
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Figure  A1 .  Finite  element  mesh  for  1  -D  column  analysis  (0.305  m  =  1  ft) 

a.  At  any  given  elevation,  the  soil  layer  possess  uniform  compressibility. 

b.  The  horizontal  strain  equals  zero. 

c.  The  vertical  settlement  at  any  given  elevation  is  uniform. 
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The  backfill  placement  analysis  of  this  1-D  column  will  serve  as  the  model  for  a 
Ko  stress-state.  The  computed  results  are  used  to  establish  that  the  value  assigned 
to  the  nominal  Poisson's  ratio  in  the  SOILSTRUCT-ALPHA  analysis  is  consistent 
with  Ko  equal  to  0.44. 
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Figure  A1  shows  the  SOILSTRUCT- ALPHA  finite  element  model  of  the  1-D 
soil  column  used  in  this  appendix.  The  moist  soil  column  is  55  ft  (16.76  m)  high. 
Water  pressures  are  assumed  equal  to  zero  at  the  time  of  backfilling  to  elevation 
425. 1  These  two  assumptions  are  consistent  with  the  initial  load  cases  of  the 
backfill  placement  analysis  of  the  new  RCC  McAlpine  lock  wall,  described  in  the 
main  body  of  the  report.  The  finite  element  mesh  for  the  soil  column  comprises 
22  soil  layers.  They  range  in  thickness  from  1.5  ft  (0.46  m)  thick  at  the  base  to  5 
ft  (1.5  m)  thick  at  the  top  of  the  column.  The  elevations  of  the  nodes  defining  the 
soil  elements  in  this  figure  were  the  same  elevations  as  the  nodes  defining  the 
new  backfill  in  Figure  6  (see  main  text).  Interface  elements  of  approximately  zero 
shear  stiffness  and  very  large  normal  stiffness  are  included  along  the  vertical 
faces  of  the  mesh  shown  in  Figure  A1  to  take  advantage  of  the  postprocessing 
capability  within  SOILSTRUCT-ALPHA.  SOILSTRUCT-ALPHA  computes  the 
horizontal  and  vertical  resultant  forces  and  their  points  of  application  (elevation) 
along  specified  regions  of  interface  elements  using  the  normal  and  shear  stress 
data.  The  pair  of  22  interface  element  columns  delineate  a  single  region  within 
the  backfill  in  this  backfill  placement  analysis.  These  results  expedite  the 
calculation  of  the  horizontal  earth  pressure  coefficient  (equivalent  to  an  at-rest 
coefficient  IQ,  in  the  1-D  soil  column)  for  the  finite  element  analysis  of  the  soil 
column. 

Total  (moist)  unit  weight  equal  to  126  pcf  (2,018.32  kg/cm3)  is  assigned  to  the 
soil.  Material  properties  assigned  to  the  hyperbolic  stress-strain  relationship  for 
the  soil  elements  are  given  in  Table  A1 .  The  values  for  the  parameters  listed  in 
this  table  are  typical  of  dense,  granular  backfill.  The  material  properties  assigned 
to  the  interface  elements  are  given  in  Table  A2. 

The  backfill  placement  analysis  of  the  Figure  A1  soil  column  is  conducted  in 
22  lifts  or,  equivalently,  22  load  increments  using  SOILSTRUCT-ALPHA.  The 
distributions  of  horizontal  effective  stress  and  vertical  effective  stress  with 
elevation  computed  within  the  soil  elements  after  placement  of  the  final  (22nd) 
lift  are  not  shown.  The  shear  stress  Txy  equals  zero  in  all  soil  elements.  The 
integral  over  the  height  of  the  soil  column  of  the  horizontal  effective  stresses  o'x , 
is  equal  to  the  horizontal  effective  force  Fx.  The  horizontal  effective  force  Fx, 
computed  by  SOILSTRUCT-ALPHA  using  the  normal  effective  stresses  within 
the  interface  elements,  equals  84,038  lb  per  ft  of  width  (1 ,226.4  kN  per  m).  The 
value  of  the  horizontal  earth  pressure  coefficient  Kj,  is  computed  using  Equation  1 
in  the  main  text. 


1  All  elevations  (el)  cited  herein  are  in  feet  referenced  to  the  National  Geodetic  Vertical  Datum. 
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Tab|e  A2 

Materia!  Properties  for  Interface  Elements  Comprising 
SOILSTRUCT-ALPHA  Finite  Element  Model  of  1-D  Soil  Column 


Analysis  (1  MPa/m  =  6,365.9  psf/ft) 


Material  Reqion 

■  ’/ 

k*,  psf/ft 

k.,  psf/ft 

Interface 

1.0x10* 

1 

- 1 

Note: 

Equations  for  Interface  Model 

The  normal  stress  at  the  center  of  the  interface  element  is  given  by 

=  kn  \ 

where  A„  is  the  average  relative  displacement  normal  to  the  interface  element.  The  shear  stress  at 
the  center  of  the  interface  element  is  given  by 

L  ■  ks 

where  A,  is  the  average  relative  shear  displacement  along  the  interface  element. 


The  denominator,  designated  as  the  Effective  Overburden  in  this  report,  is  the 
integral  of  the  vertical  effective  stress  distribution  for  the  55-ft-  (16.76  m)  high 
soil  column  (not  shown).  With  the  Effective  Overburden  computed  using 
Equation  3  (see  main  text)  equal  to  190,575  lb  per  ft  run  of  wall  (2,781.2  kN  per 
m),  Kj,  is  equal  to  0.44.  With  the  lateral  strain  e„  equal  to  zero  along  the  soil 
column,  K,,  (Equation  2  in  main  text)  is  equivalent  to  K^,. 

One  of  the  material  parameters  assigned  to  the  soil  is  the  value  for  the  nominal 
Poisson's  ratio.  This  nominal  Poisson's  ratio  used  in  SOILSTRUCT-ALPHA 
differs  from  the  traditional  strength  of  materials  definition  of  Poisson's  ratio  u. 
The  complete  derivation  of  and  its  corresponding  value  of  u  is  given  in 
Appendix  C  of  Ebeling  and  Pace  (1997).  Using  the  relationship 


l-d-2vj 

(1  -  sincj>)(l  -  K0)  1 

F  2  K0  sin<}> 

l+d-2vj 

(1  -  sin<|> )( 1  -  K0)  2 

F  2  K0  sint}) 

is  computed  equal  to  0.025  for  the  Table  A1  material  properties  with 
equal  to  0.436  or  0.44.  The  value  for  K,,  from  the  results  of  the  SOILSTRUCT- 
ALPHA  analysis  of  the  soil  column  are  consistent  with  the  results  from  using 
Equation  Al. 
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with  RF  =  0.7,  SL  =  0.5 1,  and  i)nom  =  0.025  results  in  u  equal  to  0.30.  SL  is 
computed  using  SOILSTRUCT- ALPHA  to  be  0.5 1 ,  on  average,  for  the  22  soil 
elements  of  Figure  Al. 

In  1978,  a  bulk  modulus  formulation  was  developed  by  Duncan  and  his 
colleagues  for  use  in  SOILSTRUCT  (Ebeling,  Peters,  and  Clough  1992).  This 
formulation  is  given  in  Table  B 1  in  this  report.  The  bulk  modulus  is  intended  to 
replace  Poisson's  ratio  as  the  second  elastic  parameter.  Calculation  of 
corresponding  (and  equivalent)  value  of  the  bulk  modulus  number  K,,  for  the  bulk 
modulus  formulation  is  made  using  Equation  CIO  in  Appendix  C  of  Ebeling  and 
Pace  (1997).  With  K  =  500,  uoom  =  0.025,  and  the  bulk  modulus  exponent  m  =  n  = 
0.5,  Kp  equals  175. 
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Appendix  B 

Comparison  of  SOILSTRUCT- 
ALPHA  Results  After 
Backfilling  to  Elevation  425  and 
Partial  Submergence  of  Backfill 
to  Elevation  395  for  Two 
Interface  Shear  Stiffness 
Models 


During  postconstruction,  partial  submergence  of  the  backfill  to  elevation 
(el)  3951  (identified  as  load  cases  25  through  38  in  Table  1  (see  main  text)),  the 
buoyancy  forces  act  upward,  unloading  the  backfill.  This  results  in  a  reversal  in 
the  direction  of  the  applied  shear  stress  increment  or,  equivalently,  unloading  of 
both  the  new  roller-compacted  concrete  (RCC)  lock-to-backfill  interface  elements 
and  the  existing  lock  concrete-to-backfill  interface  elements.  Interface  tests 
reported  in  Peterson  et  al.  (1976)2  show  that  the  interface  elements  follow  a 
steeper  shear-stress-versus-relative-shear-displacement  curve  than  that  described 
by  the  “primary”  hyperbolic  loading  curve  (the  tangent  shear  stiffness  k^  =  k^  [1  - 
RfiSLJ2 ).  To  model  this  “stiffer”  behavior,  the  interface  model  parameter  Rfi  is 
set  equal  to  zero  for  these  interface  elements  during  unloading  for  the  results 
reported  in  the  main  body  of  the  report.  Figure  Bla  illustrates  the  concept  behind 
the  stiff  unload  interface  model.  This  appendix  compares  these  results  with  those 
in  which  the  “primary”  hyperbolic  loading  curve  is  followed  during  partial 
submergence  of  the  backfill  to  el  395.  That  is,  the  “shear  softening”  is 
maintained  during  unloading  (i.e.,  the  term  Rfi  is  maintained  at  its  Table  5  (in 
main  report)  value  for  all  interface  regions).  This  soft  unload  interface  model  is 
illustrated  in  Figure  Bib. 


1  All  elevations  (el)  cited  herein  axe  in  feet  referenced  to  the  National  Geodetic  Vertical  Datum. 

2  References  cited  in  this  appendix  are  listed  in  the  References  at  the  end  of  main  text. 


Appendix  B  Comparison  of  SOILSTRUCT-ALPHA  Results 


B1 


Figure  B1 .  Two  models  for  interface  shear-stress-versus-relative-shear-displacement  behavior 


Figures  B2  and  B3  show  the  resulting  effective  normal  and  shear  stress 
distributions  along  the  RCC  lock-to-limestone  interface  after  backfilling  to  el  425 
and  partial  submergence  of  the  backfill  to  el  395  (load  case  38)  for  the  “stiff’  and 
“soft”  interface  shear-stiffness-versus-relative-shear-displacement  models  used 
during  unloading  (load  cases  25  through  38).  Both  effective  normal  and  shear 
stress  distributions  are  nonlinear.  As  anticipated,  the  greatest  values  of  stress  are 
computed  below  the  toe  of  the  new  RCC  lock  (x  =  0  ft  (m))  because  the  greatest 
concrete  mass  is  concentrated  above  this  region  and  because  of  the  overturning 
moment  about  the  toe  resulting  from  the  55  ft  (16.76  m)  of  backfill.  The  effective 
normal  stress  below  the  heel  of  the  new  RCC  lock  (x  =  47  ft  (14.33  m))  is  equal 
to  0  psf  (kPa)  and  307  psf  (14.7  kPa)  for  the  stiff  and  soft  interface  models  used 
during  unloading,  respectively.  Full  base  area  contact  is  maintained  in  both 
analyses. 

Table  B 1  summarizes  the  results  of  calculations  made  to  determine  the 
resultant  horizontal  effective  force  Fx  and  horizontal  earth  pressure  coefficient  Kj, 
at  section  A- A.  The  resultant  horizontal  effective  force  of  the  Figure  14  (see  main 
text)  distribution  of  o'x,  Fx,  is  equal  to  79,102  lb  per  ft  run  of  wall  (1,154.4  kN 
per  m)  and  acts  at  el  390.99  (=  0.38Htackfill)  along  section  A-A  when  a  stiff 
interface  shear  stiffness  model  is  used.  The  corresponding  value  for  the 
horizontal  earth  pressure  coefficient  is  0.459  (Figure  17  in  main  text).  The 
resultant  horizontal  effective  force  of  the  distribution  of  o'x  (not  shown)  is  made 
for  the  analysis  using  a  soft  interface  shear  stiffness  model  during  unloading  and 
Fx,  is  computed  equal  to  77,490  lb  per  ft  run  of  wall  (1,130.9  kN  per  m)  and  acts 
at  el  391.15  (=  0.38Htackfm)  along  section  A-A.  The  corresponding  value  for  the 
horizontal  earth  pressure  coefficient  K,,  is  0.45.  The  results  from  these  two 
SOILSTRUCT-ALPHA  analyses  show  that  effect  of  the  type  of  interface  shear 
stiffness  model  used  during  partial  submergence  or,  equivalently,  unloading,  of 
the  backfill  to  be  minor  in  terms  of  the  computed  values  for  Fx  and  K*  along 
section  A-A. 

Similarly,  Table  B 1  summarizes  the  results  of  calculations  made  to  determine 
the  resultant  vertical  shear  force  Fv  and  vertical  earth  pressure  coefficient  at 
section  A-A.  The  resultant  vertical  shear  force  of  the  Figure  15  (see  main  text) 
distribution  of  xxy ,  Fv ,  is  equal  to  6,775  lb  per  ft  run  of  wall  (98.9  kN  per  m) 
along  section  A-A  when  a  stiff  interface  shear  stiffness  model  is  used.  The 
corresponding  value  for  the  vertical  earth  pressure  coefficient  is  0.039 
(Figure  18  in  main  text).  The  resultant  vertical  shear  of  the  distribution  of  Txy  (not 
shown)  is  made  for  the  analysis  using  a  soft  interface  shear  stiffness  model  during 
unloading,  and  Fv  is  computed  equal  to  8,900  lb  per  ft  run  of  wall  (129.9  kN  per 
m)  along  section  A-A.  The  corresponding  value  for  the  horizontal  earth  pressure 
coefficient  K*  is  0.052.  The  results  from  these  two  SOILSTRUCT-ALPHA 
analyses  show  that  the  effect  of  the  type  of  interface  shear  stiffness  model  used 
during  partial  submergence  or,  equivalently,  unloading,  of  the  backfill  impacts 
the  computed  values  for  Fv  and  1C,,  along  section  A-A,  as  anticipated.  The  soft 
interface  shear  stiffness  model  overpredicts  the  values  of  Fv  and  IC,  for 
section  A-A  by  31  percent.  This  error  is  unconservative  since  Fv  is  a  stabilizing 
force. 
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Effective  Normal  Stress  (psf) 


Figure  B2.  Effective  normal  stress  distributions  along  the  base  of  new  RCC  lock  wall  -  partial 
submergence  to  el  395  (0.305  m  =  1  ft,  47.88  Pa  =  1  psf) 

Table  B2  summarizes  the  resultant  forces  and  their  points  of  application  on 
the  free  body  of  the  new  RCC  lock  wall  through  which  imaginary  section(s) 
shown  in  Figure  B4  are  made.  These  resultant  forces  and  their  points  of 
application  are  computed  by  SOILSTRUCT-ALPHA  after  backfilling  to  el  425 
and  after  partial  submergence  of  the  backfill  to  el  395  (load  case  38).  The 
Table  B2  results  show  that  the  shear  force  Fs  along  the  back  of  the  new  RCC  lock 
wall  is  increased  by  4.9  percent  with  the  use  of  the  soft  interface  shear  stiffness 
model  during  unloading  (load  cases  25  through  38).  The  downdrag  force  Fsis  a 
stabilizing  force,  acting  in  the  direction  of  the  heel  of  the  wall.  More 
importantly,  greater  effective  normal  stresses  are  computed  below  the  heel  of  the 
new  RCC  lock  when  the  soft  interface  shear  stiffness  model  is  used  during 
unloading  (refer  to  Figure  B2). 
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Figure  B3.  Shear  stress  distributions  along  base  of  new  RCC  lock  wall  -  partial 
submergence  to  el  395  (0.305  m  =  1  ft,  47.88  Pa  =  1  psf) 


Table  B1 

Summary  of  SOILSTRUCT-ALPHA  Results  Along  Section  A-A  for 
“Soft”  and  “Stiff”  Interface  Shear  Stiffnesses  During  Partial 
Submergence  of  Backfill  to  El  395  (Load  Case  38)  (14.594  N 
perm  =  1  lb  per  ft  run  of  wall) 


k*  During 
Unloading 


Effective 

Overburden,  lb  per 
ft  run  of  wall 

F„\  lb  per  ft 
run  of  wall 

Fv,  lb  per  ft 
run  of  wall 

K* 

172,325 

79,102 

6,775 

0.459 

0.039 

172,325 

77,490 

8,900 

0.450 

0.052 

Appendix  B  Comparison  of  SOILSTRUCT-ALPHA  Results 


Table  B2 

Summary  of  SOILSTRUCT-ALPHA  Results  After  Submergence  to  El  395  (Load  Case 

38)  for  "Soft"  and  "Stiff"  Interface  Shear  Stiffnesses  During  Partial  Submergence  of 
Backfill  (14.594  N  per  m  =  1  lb  per  ft  run  of  wall,  0.305  m  =  1  ft) 

Resultant  Normal  and  Shear  Forces 

k«.  During 
Unloading 

F„  lb  per  ft 
run  of  wall 

F„,  lb  per 
ft  run  of 
wall 

Ir..  ft 

WL 

&mob> 

deg 

T,  lb 

N',  lb 

X«,ft 

WB 

"Stiff"  with  R,,  =  0 

55,404 

128,436 

23.32 

0.36 

23.33 

93,820 

347,902 

16.18 

0.34 

“Soft"  Hyperbolic 
Primary  Curve 

58,098 

128,348 

23.17 

0.36 

24.35 

92,298 

350,127 

16.44 

0.35 

Resultant  Water  Pressure  Forces 

Load  Case 

U. _ lb  per 

ft  run  of 
wall 

Xubwt  ft 

B 

Ufront, 

lb  per 
ft  run 
of 
wall 

■ 

U^,  lb 
per  ft 
run  of 
wall 

lubackt  ft 

38 

66,435 

25.97 

0.55 

5,272. 

8 

4.33 

23,113 

9.88 

Note: 

where  B  =  47  ft,  L  =  65.19  ft,  and  tan  (  6^  )  = 

/ 

*  n 

The  values  of  stress  level,  SI* ,  computed  within  the  interface  elements  change 
with  the  submergence  of  the  backfill.  During  the  postconstruction,  partial 
submergence  of  the  backfill  to  el  395,  the  buoyancy  forces  act  upward  unloading 
the  backfill.  This  results  in  a  reversal  in  the  direction  of  the  applied  shear  stress 
increment  or,  equivalently,  unloading  of  the  22  new  RCC  lock-to-backfill 
interface  elements.  This  is  accompanied  by  a  reduction  in  the  effective  stresses 
normal  to  the  back  of  the  lock  wall  within  those  interface  elements  located  below 
the  water  table  (el  395).  At  the  end  of  backfilling  to  el  425  (load  case  24),  the 
average  value  of  stress  level  SLj  for  the  interface  elements  equals  0.70.  After 
partial  submergence  of  the  backfill  to  el  395  (load  case  38),  the  results  of  the 
SOILSTRUCT-ALPHA  interaction  analysis  using  the  stiffer  inter-face  model 
during  unloading  (Figure  Bla)  shows  that  the  average  reduction  in  stress  level  SL, 
for  the  interface  elements  is  8.6  percent  (ranging  from  zero  to  a  23-percent 
reduction).  Changes  in  both  shear  and  normal  effective  stresses  during 
submergence  account  for  the  reduction  in  values  of  the  stress  levels  SL,  for  the  22 
interface  elements. 
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Figure  B4.  Planes  along  surfaces  of  the  new  RCC  lock  on  which  resultant  forces  are  computed  -  submerged 
to  el  396  (0.305  m  =  1  ft) 
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A  new  lock  is  to  be  constructed  at  the  McAlpine  Locks  on  the  Ohio  River  at  Louisville,  KY .  The  new  lock  chamber  will  be 
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SOILSTRUCT- ALPHA  is  used  to  analyze  the  complex  interactions  among  the  lock,  the  backfill,  and  the  foundation  rock  strata. 

The  analysis  is  accomplished  in  three  phases.  The  first  phase  of  the  analysis  introduces  the  self-weight  of  the  existing  locks 
and  existing  backfill  to  the  shale  and  limestone  foundation.  The  second  phase  models  the  incremental  construction  of  the  RCC 
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interaction  analysis  of  the  north  wall  section  of  the  new  McAlpine  Lock  is  modeled  using  38  stages  of  loading  in  the 
SOILSTRUCT- ALPHA  analysis.  The  results  of  this  analysis  are  described  in  detail  in  this  report. 


